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ABSTRACT 

We present new observational results on the kinematics of the damped Lya systems. Our full sample 
is now comprised of 31 low-ion profiles and exhibits similar characteristics to the sample from Paper I. 
The primary exception is that the new distribution of velocity widths includes values out to a maximum 
of nearly 300 km s"^ , ~ 100 km s~^ greater than the previous maximum. These high velocity width 
systems will significantly leverage models introduced to explain the damped Lya systems. Comparing 
the characteristics from low-redshift and high-redshift sub-samples, we find no evidence for significant 
evolution in the kinematic properties of protogalaxies from z — 2.0 — 3.3. 

The new observations give greater statistical significance to the main conclusions of our first paper. 
In particular, those models inconsistent with the damped Lya observations in Paper I are ruled out at 
even higher levels of confidence. At the same time, the observations are consistent with a population of 
rapidly rotating, thick disks (the TRD model) at high redshift, as predicted by cosmologies with early 
structure formation. 

Buoyed by the success of the TRD model, we investigate it more closely by considering more realistic 
disk properties. Our goal is to demonstrate the statistical power of the damped Lya observations by 
investigating the robustness of the TRD model. In particular, we study the effects of warping, realistic 
rotation curves, and photoionization on the kinematics of disks in the TRD model. The principal results 
are: (1) disk warping has only minimal effect on the kinematic results, primarily influencing the effective 
disk thickness, (2) the TRD model is robust to more realistic rotation curves; we point out, however, 
that the rotation curve derived from centrifugal equilibrium with HI gas alone does not yield acceptable 
results, rather flat rotation curves such as those generated by dark matter halos are required, and (3) 
the effects of photoionization require thicker disks to give consistent velocity width distributions. 

Subject headings: cosmology — galaxies: evolution — galaxies: quasars — absorption lines 



1. INTRODUCTION 

This paper marks the third in a series of papers on the 
kinematics of the damped Lya protogalaxies. These HI 
gas layers observed along sightlines to distant QSO's are 
widely b elieved to be the gaseous progenitors of modern 
galaxies ( Wolfe 1995 ). Hence an examination of damped 
systems at high redshift provides insight into the pro- 
cess of galaxy formation. For instance, identifying the 
physical nature of these systems may distinguish between 
the monolithic collapse model ( Eggen, Lyndcn-BcU , & 



Sandage 1962 ) and the hierarchical scenar io favored by 
standard cos mogony. In our first paper (Prochaska & 



Wolfe 1997b, hereafter PW), we demonstrated that the 



kinematics of damped Lya systems at high redshift are 
consistent with these systems being thick, rapidly rotat- 
ing disks; it is a description not unlike that predicted 
by monolithic collapse formation scenarios. At the same 
time, we found damped Lya systems cannot be simple 
exponential disks in a cluster normalized Stand ard Cold 
Dark Matter cosmology (e.g. Kauffmann 1996| ). Subse- 
quently, Jedamzik and Prochaska (1998) tightened this 
conclusion by considering a range of disk characteristics 
and CDM normalizations. They found that only a finely 
tuned disk model within the framework of CDM could be 
made marginally consistent with the damped Lya obser- 

^Visiting Astronomer, W.M. Keck Telescope. The Keck Obser- 
vatory is a joint facility of the University of California and the Cali- 
fornia Institute of Technology. 



vations. Recently, Haehnelt et al. (1997) have offered an 
alternative model for damped systems as gaseous proto- 
galactic clumps undergoing infall within dark matter ha- 
los which may be consistent with the kinematic charac- 
teristics of the damped Lya systems. Such a description 
lends itself naturally to the hierarchical cosmologies where 
merging plays a vital role. A future paper will address this 
model in greater detail. For the present work, we will fo- 
cus on the interpretation of damped Lya systems as thick 
rotating disks at high redshift. 

In PW, we analyzed the low-ion profiles from 17 damped 
Lya systems and compared their kinematic characteristics 
with those of simulated profiles derived from several phys- 
ical models. Of the models tested, we found the thick 
rotating disk model to be the only model consistent with 
the observations. The basic assumptions of the model are 
a flat rotation curve and an exponential gas distribution, 
both chosen to roughly correspond with the observations 
of local spiral galaxies. In this paper we present low-ion 
profiles for 14 additional damped systems. Therefore, our 
full kinematic sample consists of 31 low-ion profiles. We 
interpret the sample with disk models containing more 
physically realistic characteristics, e.g. disk warping and 
photoionization. We have two primary goals in mind: (1) 
to test the robustness of the interpretation of damped Lya 
systems as disks and (2) to determine the effects on our 
conclusions regarding the thickness and rotation speed of 
these disks. 



1 



2 



PROTOGALACTIC DISK MODELS 



In §2 we review the terminology and methodology in- 
troduced in PW. The new data are presented and tested 
against the models of PW in §3. We investigate the effects 
of more realistic disk properties in §4 and in §5 we present 
a summary. 

2. A REVIEW 

Our strategy is to compare the kinematic characteristics 
of low-ion profiles in damped Lya systems with the simu- 
lated profiles derived from various physical models. In this 
manner, then, we investigate the agreement between the 
models and the damped systems. In this section we review 
the methodology and the thick rotating disk (TRD) model 
presented in PW. 

To characterize the kinematics of the damped systems 
we focus on low-ion profiles such as Sill 1808. After 
normalizing them to unit continuum strength, we cre- 
ate a binned apparent optical depth array, f{v), where 
f(t;) =< ln[/c//(u)] >, I{v) is the intensity at velocity 
V, Ic is the continuum intensity, and the average is taken 
over « two resolution elements centered at v. We focus on 
f(v) rather than I{v) to minimize the effects of saturation, 
and we smooth in optical depth to both limit the effects 
of Poisson noise and focus on features as opposed to in- 
dividual components. Each profile is characterized by its 
signal-to-noise ratio (SNR), instrumental resolution, and 
peak optical depth ratio, f{vpk)/cr{f), with Vpk the posi- 
tion of the maximum of f{v) in velocity space. 

In PW, the only model considered which is consistent 
with the damped Lya kinematic characteristics is the TRD 
model. For this model, we assume: 



An exponential gas distribution 



n(R, Z) = no exp 



R_ 

Rd 



\Z\ 
h 



(1) 



where R and Z are cylindrical radius and vertical dis- 
placement from midplane, tiq is the central gas den- 
sity, Rd is the radial scale length, and h is the vertical 
scale height. Note no is related to the perpendicular 
column density at i? = by N^iO) — 2noh. 

• A flat rotation curve parameterized by the rotation 
speed, Vrot 

• A random Gaussian velocity field parameterized by 
a 1-d velocity dispersion, CTcc- 

For the model systems considered here, we derive simu- 
lated profiles with identical SNR, instrumental resolution, 
and peak optical depth as the empirical data set. We em- 
ploy the same Monte Carlo techniques as those used in 
§3.2 of PW. However, we now consider 31 damped Lya 
systems which includes the 17 systems used in PW. To suf- 
ficiently sample the model distributions requires « 10000 
sightlines, thus we perform 31 runs of 333 sightlines for 
every simulation. The SNR, f{vpk) /(j{f), and instrumen- 
tal resolution of the simulated profiles for each set of 333 
sightlines is varied accordingly. 

We quantify the kinematic characteristics of the profiles 
with 4 test statistics: 



• Velocity Interval Test - The width of the profiles in 
velocity space, (/av) or Av. Formally, we calculate 
the total optical depth, Ttot of the profile in velocity 
space and define Av to be the velocity width encom- 
passing 0.9 Ttot by trimming 0.05 not off each edge of 
the profile. 



• Mean-Median Test 
profile 

fmi 



Measures the asymmetry of the 



I'^mean '^median] 



(2) 



Edge-Leading Test - Designates how edge-leading 
the peak of the profile is 



fedg 



\Vpk 



(/av/2) 



(3) 



• Two Peak Test - Designates how edge-leading the 
2nd peak of the profile is 



f2pk = ± 



\V2pk 



(/av/2) 



(4) 



In the above equations, Vpk and V2pk indicate the position 
of the peak and second strongest peak of the f{v) pro- 
file, and Vmean and Vmedian are the mean and median of 
the f{v) profile. We perform the statistical tests on the 
damped Lya and simulated profiles in the same manner 
and compare the resulting distributions with the two-sided 
Kolmogorov-Smirnov (KS) Test. 

3. DAMPED LYA PROFILES 

3.1. DATA SAMPLES AND SELECTION CRITERIA 

The full kinematic data sample, designated Sample B, 
comprises 17 damped Lya systems from PW (Sample A) 
and 14 new systems. Table |^ contains the journal of 
the 14 new observations. In Table |^ we list the follow- 
ing properties of the 14 new profiles: absorption redshift 
(column 2), log of iV(III) (column 3), Fe abundance, i.e., 
log(Fe/II) — log(Fe/H)o (column 4), Zn abundance (col- 
umn 5), the low- ion transition (column 6), and f{vpk) / (j{f) 
(column 7). From each new damped Lya system we se- 
lect a low-ion transition satisfying strict selection criteria. 
However, in order to include the highest redshift systems, 
we relax some of the selection criteria established in PW. 
This is because (1) the data have poorer SNR and (2) the 
high z systems have lower metallicity. Thus, even strong 
transitions (e.g. Fe II 1608) are difficult to measure; for 
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Table 1 
QSO AND Observational Data 



QSO Date Exp z^m Res SNR Data 







(s) 




(km/s) 




set 


Q0019-01 


F96 


35000 


4.528 


7.5 


18 




Q0347-38 


F96 


12600 


3.23 


7.5 


33 


W 


Q0951-04 


S97 


30600 


4.369 


7.5 


13 


W 


Q1005+36 


S97 


3600 


3.17 


7.5 


12 




Q1346-03 


S97 


31000 


3.992 


7.5 


29 


w 


Q1759+40 


F96 


10400 


3.05 


6.6 


33 


w 


Q2348-14 


F96 


9000 


2.940 


7.5 


41 


w 


Q0930+28 


S97 


17000 


3.436 


6.6 


26 


SL^ 


Q1104-18 


S97 


19500 


2.32 


6.6 


40 


SL 


Q1850+40 


S97 


11880 


2.120 


6.6 


16 


SL 


Q2233+13 


S97 


15000 


3.274 


6.6 


16 


SL 


Q2343+34 


S97 


40000 


2.58 


6.6 


35 


SL 


Q2344+12 


S97 


15000 


2.763 


6.6 


29 


SL 



"W - Data acquired by Wolfe & Prochaska 
''QSO coordinates kindly provided by Becker (1998) 
^SL - Data kindly provided by W. L. W. Sargent, Limin Lu, and collaborators 



Table 2 



Additional Sample of DLA Systems 



QSO 




log[Af(HI)] 


[Fe/H] 


[Zn/H] 


Transition 


TiVpk/crr) 


Q0019-01 


3.439 


20.9 


> -1.7^^ 


N/A 


Fell 1608 


b 


Q0347-38 


3.025 


20.8 


-1.8 


N/A 


Fell 1608 


33 


Q0951-04A 


3.859 


20.6 


-2.04 


N/A 


Sill 1526 


b 


Q0951-04B 


4.203 


20.4 


-2.94 


N/A 


Sill 1190 


10 


Q1005+36 


2.799 


20.6 


-2.0 


N/A 


Fell 1608 


13 


Q1346-03 


3.736 


20.7 


< -1.7 


N/A 


Sill 1304 


29 


Q1759+40 


2.625 


20.8 


-1.27 


N/A 


Sill 1808 


41 


Q2348-14 


2.279 


20.6 


-2.35 


< -1.19 


Fell 1608 


40 


Q0930+28 


3.235 


20.3 


-2.03 


N/A 


Fell 1608 


31 


Q1104-18 


1.661 


20.8 


-1.40 


-0.865 


Sill 1808 


44 


Q1850+40 


1.990 


21.4 


-1.33 


-0.70 


Znll 2026 


16 


Q2233+13 


3.149 


20.0 


-1.5 


N/A 


Fell 1608 


19 


Q2343+34 


2.430 


20.3 


-1.20 


N/A 


Sill 1808 


25 


Q2344+12 


2.538 


20.4 


-2.0 


N/A 


AlII 1670 


36 



"^AU limits are 3cr limits 
Profile is saturated at the peak 
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Fig. 1. — Velocity profiles of low- ion transitions for the 14 new damped Lya systems. For each profile, v = km s ^ corresponds to the 
redshift labeled in the plot. 
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Fig. 2. — Average optical depth profiles (binned over 9 pixels) for the 31 transitions comprising our new total empirical data set. The solid 
vertical lines designate the velocity interval, Av, while the dash-dot lines indicate peaks and the dotted line is the 2cr(f ) array. The profiles 
have been shifted in velocity space such that the left edge of the velocity interval coincides with i; = and reflected in about half of the cases 
so that the strongest peak always lies on the left edge of the interval. 
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several of these damped systems no profile satisfies all 
the selection criteria of PW. Because we wish to max- 
imize f(v)/a(f) at the edge of each profile, in partic- 
ular to ensure an accurate measurement of the velocity 
width, we adopt several profiles with peak normalized flux 
I{vpk)IIc < 0.1, which violates the previously established 
selection criterion (i.e. as in PW). In addition, we have re- 
laxed the criterion that f{vp)^) / a{f) > 20 on the grounds 
that numerical tests indicate this criterion was too conser- 
vative. In Figure |l| we present the velocity profiles as nor- 
malized intensity versus velocity, and the corresponding 
binned optical depth arrays in Figure |[ Note that many 
of the new profiles (e.g. Q0019-01, Q1104-18) exhibit 
the same edge-leading asymmetry characteristic observed 
for the profiles of Sample A. Furthermore, the distribu- 
tion of profile widths resembles that of Sample A, but the 
new distribution extends to 290 km s~^ which exceeds the 
200 km s~^ maximum width in Sample A. In fact, three 
of the 14 profiles have Av > 200 km s~^. The incidence 
of Av > 200 km s~^ in only one of 17 previous velocity 
profiles is due to small number statistics. 



(a) qool 



(b) Q0951A 




Relative Velocity (km/s) 



Fig. 3. — Comparison of a mildly saturated low-ion profile with 
a weak profile for the damped hya systems toward (a) Q0019— 01 
at z = 3.439 and (b) Q0951-04 at z = 3.859. Note the profiles 
track one another very closely and therefore exhibit nearly the same 
kinematic characteristics, irrespective of the saturation. 



Sample B includes two profiles with I{v)/Ic < 0.1, in- 
cluding one with I{v)/Ic < 0.0 where the flux falls below 
zero due to the effects of sky subtraction and Poisson noise. 
As emphasized in PW, the statistical tests used to charac- 
terize the profiles are sensitive to saturation; for example, 
a false determination of Vpk leads to incorrect values for 
all 4 statistical tests. We have strong reason to believe, 
however, that we are not contaminating the results by in- 
cluding these two profiles. First, in each case the profile 
has I(y)/Ic < 0.1 for only a few pixels. Therefore, the cen- 
troid of the peak is accurately determined to within a few 
km s^^ . Secondly, we can compare the saturated profiles 
with weaker, noisier, profiles from the same damped sys- 
tem. In Figure 0, we show this comparison. Note that the 



profiles track one another very closely exhibiting nearly 
identical kinematic characteristics. Therefore, the kine- 
matic results are not compromised by the minimal satu- 
ration evidenced in these two profiles. To calculate f{v) 
over the interval where I{v)/Ic < 0.0 we arbitrarily set 
I{v)/Ic = o'(^)/2 and calculate f{v) accordingly. Having 
performed numerical experiments which demonstrate that 
the original criteria were too strict, we establish a new set 
of saturation criteria: (1) the line profile must be free of 
blending with other absorption line profiles, (2) the profile 
must not saturate for an excess of 1 resolution element, 
(3) the saturated region must not exceed one-fourth of the 
profile velocity width, (4) f{vpk) / cr{f) > 10 instead of the 
previous value of 20. Together these criteria prevent an 
inaccurate determination of Vpk and any resulting error in 
the test statistics. 

To facilitate any parallel analysis performed with the 
damped Lya su rveys, where an ac curate iV(III) is always 
measured (e.g. Wolfe et al. 1995), and to allow for the 



most accurate comparison with the Monte Carlo simula- 
tions, we construct a subset of Sample B by imposing a 
stricter (HI) criterion than that of PW. We refer to this 
subset of Sample B as Sample C. The criterion requires 
an A^(HI) measurement for every damped system which 
must exceed N thresh = 2 x 10^" cm^^. This eliminates 
two systems from Sample A (Q1946+60A and Q0449-13) 
inferred to be damped Lya systems on account of their 
large ionic column densities and one system (Q2212— 16) 
with a measured Af(HI) < Nthresh- In addition, the crite- 
rion eliminates Q2233+13 from the new systems. There- 
fore, Sample C contains profiles from 27 damped systems. 
To investigate the evolution of the kinematic characteris- 
tics of the damped Lya systems with redshift, we divide 
Sample C at its median redshift and consider two sub- 
sets: (i) a low redshift set of profiles with z = 2.06 and 
Zmedian = 1.96 referred to as Sample D and (ii) a high 
redshift set of profiles with z — 3.24 and Zmedian — 3.35 
referred to as Sample E. Tabic ^ summarizes the 5 data 
samples. 

Table 3 
Data Samples 



Sample 


Nsys 


Comment 


A 


17 


Original selection criteria 


B 


31 


New selection criteria 


C 


27 


Af(HI) measurement required 


D 


14 


Low Redshift Cut of Sample C 


E 


13 


High Redshift Cut of Sample C 



Figure ^ shows the statistical test distributions for Sam- 
ples A— E. The results for Samples A and B are in good 
agreement indicating that the addition of the new systems 
will serve to strengthen the conclusions of PW. The only 
significant difference between the two samples is the tail 
observed out to nearly 300 km s~^ in the /av distribu- 
tion of Sample B. One also notes Samples B and C exhibit 
nearly identical distributions, therefore the stricter A^(HI) 
criterion has little effect on our results. Similarly, the only 
disagreement observed by eye between Samples D and E 
is in the Mean-Median Test, but statistically the two dis- 
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tributions are consistent {Pks = 0.6). We argue, there- 
fore, that there is httle evidence for an evolution of the 
kinematic characteristics with redshift. This will have sig- 
nificant impact on interpreting the damped Lya systems 
in terms of the competing galaxy formation scenarios and 
will be further addressed in a future paper. 



(a) Velocity Test 



(b) Mean-Median Test 



Saniple B 



Sample A 

1,1, u 



Sample C 



1^ 



Sample D 

PI J 



Sample E 



4 r 

3 

z 
1 



(c) Edge -Leading Test 



: 1 1 1 1 

[ Sample A 

' P inpi 


■ 1 1 

XX- 




[ Sample B 






[ Sample C 






[ Sample D 

M |l , i| 






[ Sample E 

r , li III 


■ , 1 r 





Sample A 



Sample B 



Sample C 



i,L ii 11 



Sample D 



III I 11 \ 1,1 |i II 



Sample E 



(d) Two-Peak Test 



Sample A 



II, I, flu |l,l, 



Sample B 



I n i r lJpll 



i 



Sample C 



I IPI \\ 



Sample D 



4^ 



I, I M III |il 1 1 



Sample E 



III 



Fig. 4. — Statistical test distributions for data Sample (A) the 
17 damped Lyo profiles from PW, (B) the 31 profiles comprising 
the current data sample, (C) 27 of the current profiles which satisfy 
stricter Af (HI) criteria, (D) a low redshift sub-sample and (E) a high 
redshift sub-sample. 



3.2. ORIGINAL MODELS REVISITED 

Figure |^ presents the statistical test distributions for 
Sample C as well as the distributions for 7 of the Monte 
Carlo models presented in PW. The Pks values are the 
KS Test probabilities that the model distribution and the 
corresponding distribution from Sample C could have been 
drawn from the same parent population. The TRDl model 
is the 'best fit' model to Sample A. For this model to 
h/Rd = 0.3, iV_L(0) = 1021-2 cm-2 and Vrot = 225 km s-\ 
where h is the vertical scale length, Rd is the radial scale 
length, A''^ (0) is the central column density normal to the 
disk, and Vrot is the fiat rotation speed. The TRD2 model 
is identical to the TRDl model except Vrot — 300 km . 
The remaining models are identical to those from §6 of 
PW. 



(a) Velocity Test 



(b) Mean-Median Test 




(c) Edge— Leading Test 



(d) Two-Peak Test 
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0,010 



0.005 
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1 1- 
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400 

1200 
800 
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Pks < OOOl 
till 



SPH 
Pug < 0.001 

I "e^ 

P„ < 0.001 



0.6 



0.8 



Fig. 5. — The statistical test distributions for the models in- 
troduced in PW compared against Sample C. Note that all of the 
models which were ruled out in PW (Pks < 5%) are ruled out at 
even higher confidence levels now. 



The addition of the new profiles has tightened the pri- 
mary results of PW in every case. In particular, we 
point out the continued failure of every model except 
the TRD model to reproduce the damped Lya obser- 
vations. The TRDl model is generally consistent with 
Sample C, but because the model adopts the same ro- 
tation speed [vrot = 225 km s^^) for all of the disks, 
it predicts no damped Lya systems with velocity width 
Aw > 250 kms~i. Hence, this simple model cannot re- 
produce the high Av tail of the new empirical /av dis- 
tribution. By contrast, the distribution from the TRD2 
model extends to Av ~ 300 km s^^ and therefore provides 
a better match to the high Av tail. Figure || plots the 
relative likelihood ratio test results for Sample C against 
the TRD model with the central HI column density (a) 
N{m) = 102" «cm-2, (b) N{m) = 1021-2 cm-2, and (c) 
A^(HI) = 102i-®cm^2 The figure indicates a lower limit 
of Vrot > 250 km s^i and that the optimal rotation speed 
is Vrot ~ 300 km s^i. However, if damped systems evolve 
to present galaxies, then a single population of disks with 
Vrot = 300 km s~i is unacceptable as such large rotation 
speeds are rarely observed. One concludes the single popu- 
lation disk model is too simple and must be revised. Mod- 
els including distributions of rotation speeds and sightline 
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encounters with multiple disks will be considered in future 
papers. In this paper we retain the assumption of a sin- 
gle population of disks and study how more realistic disk 
properties affect the test statistic distributions. 




(a) 



(b) 





400 r 




300 - 






s 


200 - 








100 - 



Fig. 6. — Isoprobability contours derived from the Likelihood 
Ratio Test for the TRD model against the f^v distribution from 
Sample C with N^{0) = (a) lO^O-^cm^, (b) lO^l ^cm^, and (c) 
1021-6 cm2. Contour levels are drawn at P = 0.01, 0.05 and 0.32. 
Note that the f/^v distribution would imply rotation speeds v > 
250 km s~l given the assumption of a single population of disks. 



4. IMPROVED DISK MODELS 

In this section we improve on the TRD model by in- 
troducing more realistic disk characteristics. We test the 
TRD model against the damped Lya observations in the 
light of these properties and thereby gauge the robustness 
of the model. The discussion focuses primarily on the ve- 
locity width test statistic because the other test statistics 
are less sensitive to changes in disk properties. 

4.1. THE WARPED DISK MODEL 

HI 21 cm observations of local spiral galaxies (e.g. 
Briggs 1990 ) suggest a significant fraction of disks are 
warped. Therefore, if the damped Lya systems evolve into 
disk galaxies, as we have argued, the effects of warps must 
be considered. Warps produce two opposite kinematic ef- 
fects: (i) some warped disks will have a large cross-section 
to sightlines which satisfy Af(HI) > 2 x 10^° cm^^ with 



large impact parameter and therefore small Aw. For in- 
stance, a tangent to the outer edge of a warped disk may 
have very large impact parameter (small Aw) yet still sat- 
isfy the A^(HI) criterion (see Figure ^a below), and (ii) 
those sightlines which doubly penetrate the disk (at a 
warped edge in addition to the unwarped inner disk, see 
Figure ||b) will tend to yield larger Aw than if they pen- 
etrated an unwarped disk. If the latter effect dominates, 
one expects a warped disk to yield a /av distribution with 
a greater number of large Aw values, which imitates a 
thicker and/or more rapidly rotating unwarped disk. In 
some cases, however, we find that the former effect dom- 
inates the results thereby lowering the median of the /av 
distribution. 

Unlike the TRD model, the warped disk model cannot 
be treated analytically. Standard modeling techniques in- 
volve decomposing the disk into a number of concentric 
rings, each with a unique orientation, with the configura- 
tion described by a 'line of nodes'. Figure ^ is a polar 
plot of a 'line of nodes' where each successive point identi- 
fies rings with larger radii. We have used different symbols 
to help distinguish the points. The radial distance (a) to 
each symbol indicates the tilt angle, and the azimuthal 
position (/3) designates the position angle of the tilt axis. 
We define the reference plane by the orientation of the 
rings with smallest radius, i.e. the inner disk. Figures |7|d- 
d show 3 views of the warped disk corresponding to this 
line of nodes. In the following, we consider 4 lines of nodes 
characterized by Rh, an effective Holmberg Radius, which 
sets the radial separation of the rings. Figure ^ plots the 
lines of nodes with respect to Rh- Included is a NULL 
line of nodes (Figure ^) describing an unwarped disk, the 
lines of node s for M83 and M33 (Figure ||b,c respectively; 
Briggs 1990| ) and an artificially constructed line of nodes 



(CLON). Because we treat the warp model numerically, 
we will compare against the results of the null models and 
not the analytic TRD Model in order to properly measure 
the consequences of warping without being biased by any 
numerical effects. Figure ^ illustrates the two kinematic 
effects mentioned above. In Figure ^ a sightline pene- 
trates a disk described by the WRP4 model and has a 
large impact parameter. The column density derived for 
this sightline is A^(HI) = 1.6 x 10^^ cm"-^ and the result- 
ing profile is plotted in the lower right panel. Note that 
this system would have A^(HI) < 2 x 10^°cm~^ for all of 
the other lines of nodes considered here and would there- 
fore not contribute to the test statistic distributions. The 
more common kinematic effect is illustrated in Figure ^ 
where we plot a sightline doubly penetrating a disk from 
the WRP2 model. The resulting profile is significantly 
wider than if the disk were not warped. 

We performed a series of Monte Carlo simulations for 
warped disks with a variety of thicknesses, Ru/Rd values, 
and lines of nodes. Figure nfl plots the /av distribution for 
8 representative warp models corresponding to the param- 
eters listed in Table ^. All of the models assume a rotation 
speed of Vrot = 250 km s~^ and a central column density, 
A^_i(0) = 10^^'^ cm~^. In general, warping yields a greater 
number of moderate Aw values (Aw w jVrot) and a few 
more large Aw values. This is the case for the M33 and 
CLON hues of nodes (e.g. compare the null model, WRPl, 
with WRP2 and WRP3) where the results resemble those 
for a thicker disk. This is not true for all lines of nodes. 
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Fig. 7. — The polar plot depicts the line of nodes for a warped disk. Each successive point on the solid line indicates rings of monotonically 
increasing radii. The azimuthal angle (/3) indicates the rotation angle of the ring and the 'radial' angle (a) designates the tilt angle of the 
ring. The remaining plots show the warped disk as a series of rings viewed from three perspectives. 
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Fig. 9. — Illustration of the two principal kinematic effects associated with warped disks: (a) a sightline with large impact parameter 
yielding a small At; which would not be included in a treatment of an unwarped disk as Af(HI) would be less than Nti^^^gf^; (b) a sightline 
doubly penetrating the disk resulting in a significantly larger Av than if the disk were not warped. The + signs indicate the position of clouds 
along the sightline. 



PROCHASKA & WOLFE 



11 



Table 4 



Warp Parameters 



Label 


h/Rd 


LON 


Rn/Rd 


WRPl 


0.2 


NULL 




WRP2 


0.2 


M33 


2.0 


WRP3 


0.2 


CLON 


2.0 


WRP4 


0.2 


M83 


2.0 


WRP5 


0.2 


CLON 


3.0 


WRP6 


0.3 


NULL 




WRP7 


0.3 


M83 


2.0 


WRP8 


0.3 


CLON 


2.0 



Note the difference between the resuhs for models WRP3 
and WRP4 which correspond to the CLON and M83 Hnes 
of nodes respectively. While both warps attain a tilt angle 
of « 40° at Rh — '^Rd, the M83 warp has a much larger 
cross-section to sightlines satisfying the Af(HI) criterion. 
As demonstrated in PW, those sightlines with large impact 
parameters tend to have smaller Av, therefore the WRP4 
model actually does more poorly than the unwarped disk, 
even though there are a number of sightlines which dou- 




bly penetrate the disk. The results are also sensitive to the 
adopted value of Rh / Rd- Comparing models WRP3 and 
WRP5 one notes the reduced effects of warping for larger 
Rh / Rd- As typical values of Rn/Rd in local disks gener- 
ally exceed 3, we expect warping to have even less of an 
effect on the /av distribution than presented here. Curves 
WRP6— 8 demonstrate the results for warping for thicker 
disks. Qualitatively, the changes in the /av distribution 
are the same as those for the thinner disks. 

To summarize, warping has only a moderate effect on 
the profile kinematics. For those lines of nodes where the 
dominant effect is of sightlines that penetrate the disk two 
times, warping mimics thicker disks, i.e. one observes a 
larger number of moderate Aw values in the /av distri- 
bution. Furthermore, warping can also lead to signifi- 
cantly fewer moderate Aw values when the cross-section 
to damped Lya system extends to higher impact param- 
eter. This is the case for systems where the outer rings 
of the warp have /3 « 180°. Having performed simula- 
tions for a large range of h, Rh /Rd, and lines of nodes, we 
quantify the results as follows: (1) in extreme cases, warp- 
ing mimics disks with up to 50% larger or smaller effective 
thickness {h/Rd value), (2) warping lends to few additional 
large Aw values in the /av distribution and therefore has 
little effect on the acceptable values for Vrot, and (3) we 
find very thin (h/Rd < 0.1) warped disks are inconsistent 
with the damped Lya observations. 

4.2. ROTATION CURVES 

In PW, the velocity field of the exponential disk model 
was given by w^ = Vrot,VR = vz = 0, i.e., we assumed a 
flat rotation curve independent of radius and height above 
the midplane of the disk. This is a good zeroth order rep- 
resentation for the rotation curve of spiral disks at large 
radii. But it breaks down at small radii because of the 
singularity implied for the mass density at i? = 0. In 
more realistic models the density approaches a finite cen- 
tral value at R less then some core radius, resulting in 
— at R — 0. The independence of with Z may also 
be unrealistic. In their analysis of the kinematics of the 
ionized gas associated with CIV absorption lines. Savage 
and Sembach (1995) find evidence that w^ = Vrot up to 
\Z\ « 5 kpc above midplane. On the other hand Sancisi 
(1998) finds evidence for a decrease in with increasing 
\Z\ in his HI studies of edge-on spirals. In fact, unless there 
is strong coupling between layers of adjacent Z this is what 
one expects. Therefore, a more physical rotation curve will 
(i) approach Vrot — km s^^ at /? = where the enclosed 
mass density presumably approaches a finite value and (ii) 
likely decrease in speed with increasing height above the 
midplane. In this section, we consider a range of rotation 
curves appropriate for a system comprised of a thick disk, 
bulge, and dark matter halo. 

4.2.1. Thick Disk 



Fig. 10. — The f^y distribution for 8 warp models compared 
against that for Sample C. The results demonstrate that warping 
has only a minor effect on the TRD Model. 



As demonstrated in PW (reaffirmed in §3.2), the kine- 
matics of the damped Lya systems are consistent with 
thick rotating disks. Furthermore, we showed that the ver- 
tical scale height, h, must be greater than one-tenth the 
radial scale length, Rd- Therefore, the thin disk approxi- 
mation is not applicable for deriving the rotation curve of 
these disks. In this section, we detail a numerical solution 
for the rotation curve of a thick disk by explicitly solving 
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Poisson's equation and combining the solution with the 
condition for centrifugal equilibrium. 

The Fourier-Bessel Transformation Method Our 

approach is an extension of the technique developed by 
Toomre (1963) for the thin disk solution and more recently 
applied by Casertano (1983) to the rotation curve for a 
thick disk at midplane. We calculate the disk potential by 
Fourier-Bessel transforming Poisson's equation, 

V2$ = 47rGp . (5) 
We adopt an exponential mass density of the form, 



p{R, Z) = po cxp{-R/Rd) exp{-\Z\/h) 



(6) 



where po, the central mass density, is related to the central 
surface density, Eg, and central HI column density, A^j_(0), 

by 



2h 



2h 



with p, the molecular weight relative to Hydrogen. This 
relation implicitly assumes that the mass of the disk is 
dominated by HI gas, i.e., we ignore any contribution from 
stars or molecular Hydrogen. 

Taking the Fourier-Bessel transformation of Poisson's 
equation, we find 



Hk,Z) 



-27rG 



+ 00 



eM-k\Z- CD p{kX)dC , (8) 



where 



Pik, = Po exp 



-ICI 



(9) 



CX3 

J Jo{kR) cxp{-R/Rd) R dR 



This integral is analytic ( Gradshteyn fc Ryzhik 1980 , 6.623 
[2]): 



p(k,C) — — ^ =- exp 

0F(i?,-2 + F)i 

Evaluating ^{k, Z), we find 



-ICI 



(10) 



^k,Z)^ 



A 



fc(i?-2 + fc2)f [hk+l 



"cxp(-Z//i) (11) 



■exp(-fcZ)] 



1 



hk-1 
where the constant A is 

A = 



[cxp{-Z/h) - exp{-kZ)] 



GttSo 
Rd 



(12) 



Finally, we take the inverse transform of 



$(i?, Z)= Jo{kR) $(fc, Z) k dk 



The thick disk rotation curve is therefore given by, 
vlUM Z) = -R^^{R,Z) 

oo 

= R I Ji{kR)^{k,Z)k^dk , 



(13) 



(14) 



which can be evaluated very accurately with standard nu- 
merical techniques. As the value for Eq derived from 
the damped Lya observations (i.e. A^j^(O)) is significantly 
lower than that observed in present galaxies, the disk ro- 
tation curve will peak at Vmax ^ 200 km (v^ cx Eoi?d) 
unless the HI gas contributes only a small fraction (< 30%) 
of the disk mass. While it is possible that stars may double 
the surface density in the inner part of the disk (see Wolfe] 



fc Prochaska 1998), we proceed under the assumption that 



the HI gas dominates. 

Thick Disk Curves Figure O shows the rotation curve 
at midplane, v^{R, 0), derived from Equation for 5 disks 
with varying thickness. Also plotted is the rotation curve 
derived with the thin disk approximation, normalized to 
have the same maximum rotation speed as the curve with 
h/Rd = 1 X 10~^. These two curves are in such close 
agreement that they are indistinguishable in Figure O, 
verifying the accuracy of the numerical solution. Note all 
of the other curves are normalized to have the same central 
column density, N±{0). Physically this normalization im- 
plies very different central mass densities for the different 
thickness disks (po oc N±{0)h~^). There is a decrease in 
with increasing h/Rd, as expected, because po is smaller 
and therefore the centrifugal force is smaller at a given 
radius R. Even for h/Rd — 0.3, however, the resulting 
rotation curve nearly traces that for the thin disk. 

The rotation curve at various heights in a disk with 
h/Rd — 0.3 is plotted in Figure |l^. Note the curves at 
Z — 1.67 h and at midplane differ by less than 10%, and 
even at Z = 6.67 h there is generally less than a 50% effect. 
We find similar results for disks of all thicknesses. Since 
nearly all of the gas arises within « 2h of the midplane, 
in general we expect the decrease in with increasing 
Z to have little effect on the kinematic results. In fact, 
it is possible the decrease in will actually increase the 
differential rotation along a given sightline. 

4.2.2. Bulge and Halo 

For the bul ge rotation curv e we assume the Hernquist 
Bulge model (Hernquist 1990) parameterized by the peak 
velocity, Vf,, and core radius, Rb 



,|Bulgc(?') = 4w^ 



rRb 



where we can relate Vb and Rb to the bulge mass, Mf, by 
the following expression: 



Mb = 2.7 X lO^M, 



Vb 



o 



200 km s- 



Rb 

lOOpc 



(16) 
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Fig. 11. — The rotation curve at midplane, Vcirc{R, 0), for five disks witli varying thickness, h/Rj. All of the curves are normalized to have 
the same central column density Nj_{0). There is only small dependence on thickness for the rotation curves. 
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Table 5 
Rotation Curves 



Crv 


"disk 


E 


Vb 


Rh 


Vh 


Rh 


Comment 




(km/s) 


(A^o/pc^) 


(km/s) 


{Rd) 


(km/s) 


(Rd) 




1 


250 












FLT 


2 


41 


18 










DLA disk 


3 


250 


1450 










Massive disk 


4 


250 


18 






250 


0.5 


DSK+HLO 1 


5 


250 


18 






250 


5.0 


DSK+HLO 2 


6 


250 


18 


180 


0.1 


250 


0.5 


DSK+HLO+BLG 1 


7 


250 


18 


180 


0.1 


250 


5.0 


DSK+HLO+BLG 2 



Meanwhile, for the halo we require a density profile 
which yields a flat rotation curve at large radii. We adopt 
the following density profile, 



p|Halo('') 



Ph 



(17) 



l + {r/Rhr ' 

where ph is the central density of the halo and Rh is the 
core radius. This gives the following rotation curve, 



I Halo 



(r) 



1 



Rh —1 

tan 

r 



r 

Rh 



(18) 



where is the halo rotation speed at large radii. 

As both the bulge and halo are spherically symmetric, 
the rotation curves in cylindrical coordinates are: 



vl\n„]frJR, Z) 



R' 



HRb 



Z"- [v'R^ + Z^ + Rb)' 



and 



V^\lii,loiR, Z) ^v^ — 



i?2 



Z2 



1 - 



(20) 



Rh 



Rh 



4.2.3. Results 

The disk rotation curve is uniquely parameterized by h 
and So, where h and all other distance scales are in units 
of the radial scale length, Rd, and So is the disk surface 
density at i? = 0. The bulge and halo models are each 
parameterized by two variables which describe the slope 
of the curve (Rb, Rh) and the maximum rotation speed 
{vbiVh) a.t Z = 0. The total rotation speed, then, is given 

by 



vl\Tot{R,Z)=vl\ii,lo{R,Z) + vl\Buigo{R,Z) (21) 
+vl\Disi.{R,Z) . 

We have performed simulations for a wide range of the six 
parameters. The main results of the analysis are presented 



in Figure ^ and the rotational parameters are given in 
Table |^. Also listed is the peak rotation speed of each 
rotation curve, Vmax- For all of the models presented here, 
we assume N±{0) = lO^^ ^ cm'^ and that h/Ra = 0.3. 

To facilitate a direct comparison with the TRD model, 
we plot the results for a flat rotation curve with Vmax = 
250 km s~^. For a given value of Vmax, the flat rotation 
curve yields a greater fraction of large velocity widths 
than any other rotation curve. Therefore, it gives the 
best agreement with the empirical distribution provided 
Vmax ^ 350 km s~^, above which the distribution has too 
many large A?; values. For model CRV2, we assume a thick 
disk with So = fimpN±{0) = ISMq/^c'^ and Rd = 10 kpc 
which implies Vmax — 41 km s~^. This model corresponds 
to A^_l(0) = 102^cm^2 as inferred from the damped Lya 
observations. The value of Rd was chosen to roughly cor- 
respond to the cross-section derived for damped Lya sys- 
tems assuming a number density approximat ely equal to 



that for spiral galaxies in the present epoch { Wolfe et al 
1995). As noted above and shown in Figure [l3|, the mass 



associated with the HI disk implies a rotation curve which 
is inconsistent with the damped Lya observations. The 
results for the CRV3 model demonstrate that a disk with 
So = 144OM0/pc2 and Rd — 5 kpc is consistent, yet this 
implies that > 95% of the disk mass is in a component 
other than HI gas. It is unlikely that So approache s even 
500Af(7i/pc2 given current estimate s of the stellar (Wolfe 
fc Prochaska 199S| ) and molecular ( pe fc Bechtold 1997] ) 
baryonic fraction. Furthermore, if Sq ~ 15OOM0/pc^ in 
every damped system, this would imply a comoving den- 
sity approximately 100 times that of visible matter in the 
present universe. On the other hand if we include a mas- 
sive dark halo {vh = 250 km s^^), the resulting /av dis- 
tribution (CRV4) is consistent with the empirical distri- 
bution provided Ru/Rd <, 1- For larger values of Rh/Rd, 
the rotation curve (CRV5) rises too slowly to yield large 
Av. Finally, we investigate the effect of including a mas- 
sive bulge in models CRV6 and CRV7 where we adopt the 
same disk and halo curves from models CRV4 and CRV5. 
We find the presence of a massive bulge {vb — 180 km s~^) 
with Rb <C Rd has little impact on the /av distribution 
as the bulge affects only the inner rotation curve. To con- 
clude, the TRD model requires a rotation curve which is 
nearly flat for Rd > 1. While this can be achieved with a 
massive disk resembling that of the Milky Way, it would 
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Fig. 13. — The left hand side plots the rotation curve at midplane for disks with varying bulge, halo, and disk components. The right hand 
side plots the /av distributions for the 8 curves compared against Sample C. 
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requires a non-gaseous baryonic component with comoving 
density two orders of magnitude greater than the density 
of visible matter in current galaxies, which is implausible. 
Assuming that the mass of the disk is dominated by the 
HI gas with iVj_(0) « 10^^ cm^^ requires the presence of a 
massive dark halo with Rh <^ Rd- 

4.3. PHOTOIONIZATION 

Over the redshift range spanning our damped Lya ob- 
servations the IGM is primarily ionized, presumably by 
the ambient UV radiation field from background quasars. 
At redshift z ^ 2.5 th e UV intensity is often d escribed as 
a modified power law ( Haardt fc Madau 1996| ) 



J{v) K. Jgi2 



^^912 



(22) 



with Jgi2 « 10~^^'^ erg/s s^^ cm^^ sr^^ Hz^^ and expo- 
nent aw —1. This background flux creates ionization 
fronts in the HI disks comprising the TRD Mod el, anal o- 



gous to those observed in present spiral galaxies ( Maloney 
1993 ). Because a complete treatment involving a solution 
of the radiative transfer equation is beyond the scope of 
this paper, we address the problem with two different ap- 
proximations. In one case, we assume a sharp HI edge to 
the disks of the TRD model. In the other, we assume the 
gas is photoionizcd if its volume density is below some crit- 
ical value. The photoionization of the disk will eliminate 
from the derived statistical sample those sightlines with 
A^(HI) K, N thresh- Since these sightlines tend to penetrate 
the disk at large impact parameters which yield small Au, 
we predict a general shift in the /av distribution to larger 
velocity widths. By analogy to the results from warped 
disk models, we expect the photoionizcd disk to mimic 
thicker disks. 

First, we consider the HI Edge Model, which incorpo- 
rates a sharp "HI edge" for the exponential disk at a col- 
umn density, Nph- Using a 2D analogy to the Stromgren 
sphere argument, we find that the column density of ion- 
ized gas above the HI gas is given by 



N, 



ph 



2(j) 



arec < n > 



(23) 



where is the flux of ionizing photons 

oo 

(i.e. (/) = / dv2'Kj^/hv)^ < n > is the average hydro- 

gen volume density, and a^ec is the Case B recombina- 
tion coefficient. Given the above value for J912, with 
< n >= 1 cm^'^ we find Nph ~ 10^*cm~^. To simu- 
late photoionization of the HI disks in the TRD model, 
we systematically remove Nph of column normal to the 
disk everywhere. This leads to a sharp HI edge at Rph — 
hi[N / 2N ph\ and effectively removes 2Nph from each 
sightline which penetrates within Rph- As J912 is poorly 
constrained observationally, we consider J912 values rang- 
ing from 10~^^ ° — 10~^^'^ er g s~^ cm~^ sr~^ Hz~^, w here 
the favored value is IQ-^^-^ ( [Haardt fc Madau 1996| ). To 
consider meaningful Nph values, we relate Nph to A^^(O), 
h, and J912 with the following equation: 



Nph = 10^^ cm 



^-2 ( h/Rd 

V 0.2 



10-21.5 



10^ 



N±iO) 



(24) 

The other approximation we make to describe the pho- 
toionization of the disk is to assume the gas is predomi- 
nantly ionized below a critical volume density, Uph- For 
the TRD model, the critical density can be defined as 



^ph 



n{R^Rph,Z^O) 



(25) 



cm 



with Rph and Nph defined as above. 
Rd ^ 5 kpc, and Nph = 10^^ - 10^° 
density cutoffs of Uph ~ [10"'' — 10~^cm~'^]. We refer to 
this photoionization model as the Critical Density model. 



Nph 
h 

With h = O.SRd, 
^ this implies 
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Fig. 14. — Test statistic distributions for the removed sightlines 
from four representative runs of the two photoionization models 
considered in the text. 



It is revealing to examine the kinematic characteristics 
of those sightlines removed from the TRD statistical sam- 
ple on account of photoionization. Figure |l^ plots the 
test statistics for the eliminated sightlines from 4 repre- 
sentative runs of each photoionization model the param- 
eters of which are given in Table ^ Included in each 
panel of the /av distribution is the percentage of sight- 
lines eliminated. Note that very few sightlines are re- 
moved in a 10000 sightline run for Nph ~ 10^^ cm~^, hence 
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the effects of photoionization are minimal for small J912, 
larger N±{0) or small h/Rd- By contrast, simulations 
with Nph = 10^" cm~^ have very many sightlines removed 
so that there is a large impact on the simulations. As 
expected, the vast majority of sightlines eliminated have 
small Aw and correspondingly small /mm and /2pk values. 
Also, the results for thin disks are essentially unaffected 
by photoionization as most sightlines yield small Aw irre- 
spective of photoionization. The results are very similar 
for the two photoionization models, although the Critical 
Density model does remove a few more large Aw. 

The results presented in Figure |lj suggest photoion- 
ization will improve the agreement between models with 
small thickness or low rotation speed (i.e. models with an 
/av distribution dominated by small Aw) and the em- 
pirical data set. However, besides eliminating sightlines 
with small iV(HI), photoionization of the edges of the 
disks tends to lower the average Aw for a given sight- 
line; because of photoionization, the HI path along a given 
sightline is reduced, which reduces the differential velocity 
along the line of sight and on average lowers Aw. The 
net result is that photoionization actually worsens the 
agreement for disks of all thickness and central column 
density, contrary to our initial expectation. Figure |l^ 
presents isoprobability Pxs(Aw) contours for disks with 
(i) 7Vi(0) = 1020-8 cm-2, (ii) 7V^(0) = lO^i-^cm-^, and 
(in) iVj.(0) = 1021-6 cm-2, for a range of h/Ra and J912 
values in the (a) HI Edge model and (b) Critical Den- 
sity model. The three contours correspond to P^s (Aw) 
values of 0.01, 0.05, and 0.32 and the x marks the high- 
est Pifs(Aw) value in the explored parameter space. The 
agreement clearly worsens for increasing J912 and decreas- 
ing h/Rd- If the UV background does have J912 ~ 10~2i^ 
then photoionization will have a significant impact on the 
results derived for in the TRD Model. In particular, the 
thickness would have to exceed h/Rd > 0.2 for disks with 
Nj_{0) = 1020-** cm-2 and 0.3 for disks with larger N±{0). 
While such large values for J912 are unlikely, they are not 
ruled out by current observations. For the J912 « 10~2i-5 
the effect on the TRD model is moderate, generally re- 
quiring less than a 30% increase in the thickness of the 
disks. 

Lastly, we examine how the photoionization models af- 



Table 6 



Photoionization Models 



Label 


Model'' 


h/Rd 


log[Af(HI)] 


-log[J9i2] 


log[Afph] 


PHT 1 


HI'' 


0.3 


21.2 


22.0 


18.7 


PHT 2 


HI 


0.3 


21.2 


20.5 


20.2 


PHT 3 


HI 


0.3 


20.8 


21.5 


19.6 


PHT 4 


HI 


0.1 


21.2 


21.0 


19.2 


PHT 5 


CD= 


0.3 


21.2 


22.0 


18.7 


PHT 6 


CD 


0.3 


21.2 


20.5 


20.2 


PHT 7 


CD 


0.3 


20.8 


21.5 


19.6 


PHT 8 


CD 


0.1 


21.2 


21.0 


19.2 
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Fig. 16. — Plots of the Af(HI) vs. Av pairs for the damped Lya systems (boxes) and four photoionization models: PHTl and PHT5 where 
J912 is small and PHT2 and PHT6 where photoionization is significant. The Pks are derived from the 2-dimensional KS Test. Note that 
photoionization clearly worsens the agreement between the model and the empirical distribution. 



feet the distribution of A^(HI), Av pairs. In PW (Fig- 
ure 13) we noted the rather poor concordance between 
the TRD model and the damped Lya observations and 
suggested photoionization may improve the agreement. 
SpecificaUy, we expected photoionization to reduce the 
number of sighthnes with small Av and small A^(III) and 
thereby improve the agreement with the damped Lya ob- 
servations. Figure 16 plots the 7V(HI), Av pairs from 
Sample C (big squares) and the A^(HI), Av pairs for 4 
of the photoionization models. The PHTl and PIIT5 
models indicate the results without significant photoion- 
ization while the PIIT2 and PIIT6 models highlight its ef- 
fects. The Pks values were derived from the 2-dimcnsional 
Kolmogorov-Smirnov test (Press 1992). For the reasons 
discussed above, photoionization worsens the agreement 
between the TRD model and damped distributions. While 
models PIIT2 and PIIT6 are inconsistent at the 99% c.L, 
it should be noted, however, that one can improve the con- 
cordance by considering a range of Vrot and N±{0) values. 
Also, as noted in Wolfe & Prochaska (1998), the presence 
of many large Av, small A^(III) pairs may indicate the 
presence of a 'hole' in the inner region of the disk as ob- 
served in many local spirals. 

5. SUMMARY AND CONCLUSIONS 

We have presented new observations on the low-ion kine- 
matics of the damped Lya systems. The full sample of 31 
profiles confirms the primary conclusions of PW, in partic- 
ular, (i) models with kinematics dominated by random or 
symmetric velocity fields are inconsistent with the damped 
Lya kinematics, (ii) the TRD model, which consists of a 



population of thick, rapidly rotating disks at high z, nat- 
urally reproduces both the observed edge-leading asym- 
metry of the empirical profiles as well as the distribution 
of velocity widths, and (iii) models incorporating centrifu- 
gally supported disks within the framework of the standard 
CDM cosmology are ruled out at high levels of confidence. 
In addition, a comparison of the kinematic properties of 
profiles of the highest redshift systems (z — 3.24) with the 
lower redshift systems {z = 2.06) reveals no significant evo- 
lution in the kinematics of the damped Lya systems. This 
last observation may place strong constraints on scenar- 
ios of galaxy formation which predict significant evolution 
over this epoch. 

Presently there are two working models which explain 
the kinematic characteristics of the damped Lya systems: 
(1) the TRD model and (2) merging protogalactic clumps 
in numerical si mulations of the stand ard Cold Dark Mat- 
ter cosmology (Haehnelt ct al. 1997). In this paper we 



have focused on the TRD model. In particular, we have 
investigated the robustness of the model to including more 
realistic disk properties, specifically disk warping, physi- 
cal rotation curves and photoionization. Given the preva- 
lence of warping in local disk galaxies, we considered its 
effects on the kinematics of the disks in the TRD model. 
We found that the results of warping are dominated by 
two competing effects. Sightlines which penetrate both 
the unwarped inner disk and the warped outer disk yield 
moderately higher Av than those simply intersecting an 
unwarped disk. At the same time, however, some warped 
disks have significantly larger cross-section to sightlines 
with large impact parameters which tend to yield small 
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Av. Having considered a number of warped disks with a 
broad range of properties, we find: (i) in extreme cases, 
warping mimics disks with up to 50% larger or smaUer ef- 
fective thickness {h/Rd value), (ii) warping leads to very 
few extra large Av values in the /av distribution and 
therefore has little consequence on the acceptable values 
for Vrot, and (iii) the lower limit to h/Rd is nearly un- 
changed as we find h must be > O.lRd for both warped 
and unwarped disks. 

In PW, we assumed a flat rotation curve, = frot, 
extending from i? = ^ oo and Z = ^ oo. In this 
paper we adopted rotation curves derived from specific 
bulge, halo, and disk components. Assuming an expo- 
nential profile is a good description of the density profile 
for the damped Lya systems, we find the rotation curves 
derived from gravity generated by the HI gas alone can- 
not reproduce the empirical /av distribution. If the rota- 
tion curve is dominated by the disk, one must introduce 
another mass component (e.g. stars, molecules) to estab- 
lish consistency. At the same time we find that the rota- 
tion curve derived from a massive halo with core radius 
Rh <, Rd also yields a /av distribution consistent with the 
observations. We believe that this latter explanation is 
more plausible. We also find the presence of the bulge to 
be largely inconsequential. 

Lastly, we studied the effects of the intergalactic pho- 
toionizing background radiation on the disk kinematics. 
We made two separate approximations to model the pho- 
toionization of the disks: (a) an HI edge model where the 
disk is photoionized at radii R > Rph with Rph set by 
the intensity of the photoionizing background and the disk 
properties and (b) a Critical Density model where all gas 
with volume density n < n^h is presumed ionized. Con- 
trary to our expectations, we find that photoionization 



tends to worsen the agreement between the TRD model 
and the damped Lya observations. I'hc effect, however, is 
not large (< 30%) for the favored value of J912 — 10~^^'^, 
but for Jgi2 = 10"^^ a substantial (> 50%) increase in the 
thickness of the disks would be required. 

In summation, then, wc find the TRD model is very ro- 
bust to tests against the damped Lya observations. The 
challenge remains, however, to consistently incorporate 
this model within a cosmological framework. While the 
clump model fits naturally within the SCDM cosmology, 
it must be demonstrated that the clump model exhibits 
similar robustness to comparisons with the damped Lya 
kinematics. While Haehnelt et al. (1997) did show that the 
clump model could explain the damped Lya observations 
from PW for a single set of parameters, a formal investiga- 
tion of the full physical parameter space with meaningful 
statistics has yet to be performed. In addition, it is not 
clear how that model will change given different cosmo- 
logical parameters, e.g. an Open Universe where merging 
plays a smaller role at z « 2.5. The model must also be 
tested against the new observations, in particular the new 
/av distribution which extends to ~ 300 km s^^. Finally, 
the fact that the numerical simulations do not reproduce 
the observed properties o f modern galaxies when evolved 
to the present universe (Navarro, Frenk & White 1995) 
suggests the model may have serious inconsistencies in the 
early universe. 

In future papers we will introduce observations of the 
high-ion transitions (e.g. CIV 1548) with the aim of further 



constraining the two working models as well as advancing 
our understanding of the ionized gas associated with the 
damped Lya systems. This gas is presumed to reside in the 
halo of these protogalaxies and therefore may give more 
direct indications of the dark matter associated with the 
damped Lya systems. We also intend to consider effects 
(e.g. multiple disks) which would improve the agreement 
between the semi - analytic models of standa rd cosmology 
dKauffmann 1996| , pVIo, Mao, fc White 199?! ). 
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